Towards the understanding of transcriptional and translational regulatory complexity

Considering the same genome within every cell, the observed phenotypic diversity can only arise from highly regulated mechanisms beyond the encoded DNA sequence. We investigated several mechanisms of protein biosynthesis and analyzed DNA methylation patterns, alternative translation sites, and genomic mutations. As chromatin states are determined by epigenetic modifications and nucleosome occupancy,we conducted a structural superimposition approach between DNA methyltransferase 1 (DNMT1) and the nucleosome, which suggests that DNA methylation is dependent on accessibility of DNMT1 to nucleosome–bound DNA. Considering translation, alternative non–AUG translation initiation was observed. We developed reliable prediction models to detect these alternative start sites in a given mRNA sequence. Our tool PreTIS provides initiation confidences for all frame–independent non–cognate and AUG starts. Despite these innate factors, specific sequence variations can additionally affect a phenotype. We conduced a genome–wide analysis with millions of mutations and found an accumulation of SNPs next to transcription starts that could relate to a gene–specific regulatory signal. We also report similar conservation of canonical and alternative translation sites, highlighting the relevance of alternative mechanisms. Finally, our tool MutaNET automates variation analysis by scoring the impact of individual mutations on cell function while also integrating a gene regulatory network.Da sich in jeder Zelle die gleiche genomische Information befindet, kann die vorliegende phänotypische Vielfalt nur durch hochregulierte Mechanismen jenseits der kodierten DNA– Sequenz erklärt werden. Wir untersuchten Mechanismen der Proteinbiosynthese und analysierten DNA–Methylierungsmuster, alternative Translation und genomische Mutationen. Da die Chromatinorganisation von epigenetischen Modifikationen und Nukleosompositionen bestimmt wird, führten wir ein strukturelles Alignment zwischen DNA–Methyltransferase 1 (DNMT1) und Nukleosom durch. Dieses lässt vermuten, dass DNA–Methylierung von einer Zugänglichkeit der DNMT1 zur nukleosomalen DNA abhängt. Hinsichtlich der Translation haben wir verlässliche Vorhersagemodelle entwickelt, um alternative Starts zu identifizieren. Anhand einer mRNA–Sequenz bestimmt unser Tool PreTIS die Initiationskonfidenzen aller alternativen nicht–AUG und AUG Starts. Auch können sich Sequenzvarianten auf den Phänotyp auswirken. In einer genomweiten Untersuchung von mehreren Millionen Mutationen fanden wir eine Anreicherung von SNPs nahe des Transkriptionsstarts,welche auf ein genspezifisches regulatorisches Signal hindeuten könnte. Außerdem beobachteten wir eine ähnliche Konservierung von kanonischen und alternativen Translationsstarts, was die Relevanz alternativer Mechanismen belegt. Auch bewertet unser Tool MutaNET mit Hilfe von Scores und eines Genregulationsnetzwerkes automatisch den Einfluss einzelner Mutationen auf die Zellfunktion

SNP and indel frequencies at transcription start sites and at canonical and alternative translation initiation sites in the human genome

Single-nucleotide polymorphisms (SNPs) are the most common form of genetic variation in humans and drive phenotypic variation. Due to evolutionary conservation, SNPs and indels (insertion and deletions) are depleted in functionally important sequence elements. Recently, population-scale sequencing efforts such as the 1000 Genomes Project and the Genome of the Netherlands Project have catalogued large numbers of sequence variants. Here, we present a systematic analysis of the polymorphisms reported by these two projects in different coding and non-coding genomic elements of the human genome (intergenic regions, CpG islands, promoters, 5' UTRs, coding exons, 3' UTRs, introns, and intragenic regions). Furthermore, we were especially interested in the distribution of SNPs and indels in direct vicinity to the transcription start site (TSS) and translation start site (CSS). Thereby, we discovered an enrichment of dinucleotides CpG and CpA and an accumulation of SNPs at base position -1 relative to the TSS that involved primarily CpG and CpA dinucleotides. Genes having a CpG dinucleotide at TSS position -1 were enriched in the functional GO terms "Phosphoprotein", "Alternative splicing", and "Protein binding". Focusing on the CSS, we compared SNP patterns in the flanking regions of canonical and alternative AUG and near-cognate start sites where we considered alternative starts previously identified by experimental ribosome profiling. We observed similar conservation patterns of canonical and alternative translation start sites, which underlines the importance of alternative translation mechanisms for cellular function

SNP and indel frequencies at transcription start sites and at canonical and alternative translation initiation sites in the human genome.

Single-nucleotide polymorphisms (SNPs) are the most common form of genetic variation in humans and drive phenotypic variation. Due to evolutionary conservation, SNPs and indels (insertion and deletions) are depleted in functionally important sequence elements. Recently, population-scale sequencing efforts such as the 1000 Genomes Project and the Genome of the Netherlands Project have catalogued large numbers of sequence variants. Here, we present a systematic analysis of the polymorphisms reported by these two projects in different coding and non-coding genomic elements of the human genome (intergenic regions, CpG islands, promoters, 5' UTRs, coding exons, 3' UTRs, introns, and intragenic regions). Furthermore, we were especially interested in the distribution of SNPs and indels in direct vicinity to the transcription start site (TSS) and translation start site (CSS). Thereby, we discovered an enrichment of dinucleotides CpG and CpA and an accumulation of SNPs at base position -1 relative to the TSS that involved primarily CpG and CpA dinucleotides. Genes having a CpG dinucleotide at TSS position -1 were enriched in the functional GO terms "Phosphoprotein", "Alternative splicing", and "Protein binding". Focusing on the CSS, we compared SNP patterns in the flanking regions of canonical and alternative AUG and near-cognate start sites where we considered alternative starts previously identified by experimental ribosome profiling. We observed similar conservation patterns of canonical and alternative translation start sites, which underlines the importance of alternative translation mechanisms for cellular function

Evaluation of the Molecular Pathogenesis of Adrenocortical Tumors by Whole-Genome Sequencing

Pathogenesis of autonomous steroid secretion and adrenocortical tumorigenesis remains partially obscure. Our aim was to identify novel genetic alterations in adrenocortical adenomas (ACA) without somatic mutations in known driver genes. Whole-genome sequencing was performed on 26 ACA/blood-derived DNA pairs without driver mutations in PRKACA, GNAS and CTNNB1 genes at previous WES (ENSAT study JCEM 2016). These included 12 cortisol-producing adenomas with Cushing syndrome (CS-CPAs), 7 with mild autonomous cortisol secretion (MACS-CPAs), and 7 endocrine-inactive ACAs (EIAs). Seven adrenocortical carcinomas (ACC) were added to the cohort. We developed a bioinformatics pipeline for a comprehensive genome analysis and to reveal differences in variant distribution. Strelka, VarScan2 and ANNOVAR software and an in-house confidence score were used for variant calling and functional annotation. Combined Annotation-Dependent-Depletion (CADD) values were used to prioritize pathogenic variants. Additional focus relied on variants in pathogenically known pathways (Wnt/β-catenin, cAMP/PKA pathway). NovoBreak algorithm was applied to discover structural variations. Two hypermutated CS-CPA samples were excluded from further analysis. Using different filters, we detected variants in driver genes not observed at WES (one p.S45P in CTNNB1 and one p.R206L in PRKACA in two different CS-CPAs). In total, we report 179,830 variations (179,598 SNVs; 232 indels) throughout all samples, being more abundant in ACC (88,954) compared to ACA (CS-CPAs: 31,821; MACS-CPAs: 35,008; EIAs: 29,963). Most alterations were in intergenic (>50%), followed by intronic and ncRNA intronic regions. A total of 32 predicted pathogenic variants were found in both coding (CADD values ≥ 15) and non-coding (CADD values ≥ 5) regions. We found 3,301 possibly damaging and recurrent variants (intergenic mutations removed) (CS-CPAs: 1,463; MACS-CPAs: 1,549; EIAs: 1,268; ACC: 1,660), mostly accumulated in intronic regions. Some of these were detected in members of the Wnt/β-catenin (CS-CPAs: 6; MACS-CPAs: 2; EIA: 1) and cAMP/PKA (CS-CPAs: 6; MACS-CPAs: 7; EIA: 4) pathways (e.g. ADCY1, ADCY2, GNA13, PDE11A). We also found a slightly higher number of structural variations in EIA (3,620) and ACC (3,486) compared to CS-CPAs (977) and MACS-CPAs (2,119). In conclusion, still unrevealed genetic alterations, especially in intronic regions, may accompany early adrenal tumorigenesis and/or autonomous cortisol secretion